NASA TM X- 72800 





This Informal documentation medium Is used to provide accel erated or 
special release of technical information to selected users. The contents 
may not meet NASA format editing and publication standards, may be re- 
vised, or may be incorporated in another publication. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
UNGLEY RESEARCH CENTER, HAMPTON, VIRGINIA 23665 


< 


I 


1 

i 


1. Hpian No. 

NASA TM X- 72800 
4. Titio and SulUttN 


2 . 


G o ^ a mw tnt Ao ci ion No. 


Flutter of Laminated Plates In Supersonic Flow 


3. Racipiant't Cataloa No. 


5. Raport Data 

Nov> 1975 

6. Rarformina Organitadon Coda 


7. AutlwrU) 


a Mormins Organiiation Raport No. 


James Wayne Sawyer 

9. Rwformifig Organiiation NamawidAddrfii 
NASA Langley Research Center 
Hampton, VA 23665 


12. Spomoring Agancy Name and Addreis 

National Aeronautics and Space Administration 
Washington, DC 20546 

IS. Suppltmantafy Notei ~ 


16. Abstract 


10. Work Unit No. 

506-17-22-01 

11. Contract or Grant No. 


13. Type of Raport and Rariod Covered 

Technical Memorandum 

14. Spcniormg Agtncy Cod* 


A solution procedure has been developed using linear small deflection theory 
for the flutter of simply supported laminated plates. For such plates, the 
bending and extensional governing equations are coupled and have cross- 
stiffness terms which do not appear in classical plate theory. An extended 
Galerkin method is used to obtain approximate solutions to the governing 
equations, and the aerodynamic pressure loading used in the analysis is that 
given by linear piston theory with flow at arbitrary cross-flow angle. A 
limited parametric study has been conducted for typical laminated composite 
plates. The calculations show that both the bending-extenslonal coupling and 
the cross-stiffness terms have a large destabilizing effect on flutter. 

Since classical plate theory does not consider bending-extenslonal coupling 
and cross stiffness terms, it usually gives Inaccurate and nonconservative 
flutter boundaries for laminated plates. 


17 


Kty Word! ISuggMttd by AuthgrWI (STAR etttgory umkrIinMl 


IS. Dittrlbutton SUMmwM 


STRUCTURAL MECHAHICS 


Panel Flutter 


Unclassified ■ 

Unlimited 


Composite Plates 
Laminated Plates 





19. SKurity OMiif. M thte riportl 

20. Sieurity Omti, (of this pogi) 

21. No. of Pages 

22. Wet* 

Unclaaalfied 

Unclaaalfied 

119 

.1:5.25 


iTIw NetioMi Technical InformMion Swvice, S|ir)n|fiaM, Virainia 221S1 
*AyaHablafraffl{ 


STIF/NASA SelMtifie and Technical Informition Facility, P.O. Bon 33, Collage Park, MD 20740 










TABLE OF CONTENTS 


CHAPTER page 

TITLE ^ 

TABLE OF CONTENTS ii 

NOTATION Iv 

LIST OF TABLES vii 

LIST OF FIGURES viii 

1. INTRODUCTION 1 

1.1 Review of Literature 2 

1.1.1 Formulation of Equations ..... 3 

1.1.2 Solutions for Symmetric Plates 

1.1.3 Solutions for Angle-ply and Cross-ply Plates . 6 

1.1.^ Approximate Solutions 8 

1.1.5 Dynamic Stability 9 

1.2 Statement of Problem ..... 9 

2. DEVELOPMENT OF GOVERNING EQUATIONS 13 

2.1 Lamina Stress-strain Equations 1^ 

2.2 Lamina and Plate Constitutive Equations .... 1^ 

2.3 Plate Governing Equations 21 

2.3.1 Equations of Motion 21 

2.3.2 Compatibility Equation 26 

2,h Boundary Conditions 28 



lii 


Oi 

01 




CHAPTER PAGE 

3. APPROXIMATE SOLUTKWS OF THE FLUTTER EQUATIONS 31 

3.1 General Laminated Plates 32 

3*2 Symmetric Laminated Plates ^1 

3.3 Angle-ply Plates 

3.^ Approximate Solutions ^9 

U. RESULTS AND DISCUSSION 52 

U.l Convergence of Results 53 

1.2 Comparison vith Literature 5l 

1.3 Symmetric Plates 56 

1.3.1 Flutter Boundaries 5T 

1.3.2 Effects of Inplane Ncrmal and Shear Loads. . 58 

1.3.3 Effect of Cross-flow 59 

1.1 Angle-ply Plates 6l 

1.1.1 Flutter Boundeuries 6l 

1.1.2 Effect of Inplane Normal and Shear Loads . . 63 

1.1.3 Effect of Cross-flow 63 

1.5 General Laminated Plates 61 

1.5.1 Flutter Boxindaries 65 

1.5.2 Effect of Plate Construction 65 

1 . 5.3 Composite Stiffened Aluminum Plates .... 6T 

5 . CONCLUSIONS AND RECOMMENDATIONS 69 

REFERENCES 7l 


I 



NOTATION 


C , C 
mn* rs 


\ l*^22 


H ,H 
mn* rs 


K„. K 

X* y 


matrix of elastic areas defined by equation (2.16) 

matrix defined by equation (2.23) 

coefficient matrix (See equation (3.18).) 

plate dimensions along x and y axis* respectively 

matrix of elastic statical moment defined by equation (2.l8) 

matrix defined by equation (2.23) 

coefficient matrix (See equation (3.18).) 

Fourier series coefficients 
coefficient matrix (See eqxiation 3.180 

matrix of elastic moments of Inertia defined by equation (2.l8) 

matrix defined by equation (2.23) 

matrix of D. . with 6»0. 

matrix of D. . with 9*0 and 6*=0 
* J 

Principal Young's moduli in 1 and 2 directions, respectively 
Airy stress function, defined by equation (2.31) 
shear modulus in 1>2 planes 
total plate thickness 

distance from plate midplane to the lamina boundary as shown 
in figure 3. 

Fourier series coefficients 
integers 

number of liters in plate » p 

V V 

nondimei\aion inplane nomeJ. loads 


iv 





V 



M 




M ,M 
x’ y 

M 

xy 

m,n 


N 


x’ 



nondimenslonal Inplane shear load ■ 

hi” 

Mach number 

total number of terms In x and y direction, respectively 
bending moment resultants as defined by equation (2.9) 
twisting moment resultant as defined by equation (2.9) 
integers 

inplane normal loads per unit length in x and y 
directions, respectively 


xy 


xy 


N' , N' 
x’ y 

N , N 
x’ y 


xy 

P 

5ij 

r,8 


u®, V® 


x,y,z 


induced inplane shearing loads 
inplane shearing loads per unit length 
Induced inplane normal loads 

inplane normal resultant as defined by equation (2.8) 
inplane shear resultant as defined by equation (2.8) 
intensity of transverse load 
constitutive coefficient matrix for specifdly 
orthotropic lamina 

constitutive coefficient matrix for generally orthotropic 
lamina 

dynamic pressure of airstream 
Integers 

midplane displacements in x and y direction, respectively 
normal displacements 
cartesian coordinates 



Vi 


6 compressibility factor) 

Y plate mass per unit area 

^xy*^yz*^xz strain components 

6 lamina thickness (See figure lU.) 

£y normal strains in x, euid y directions, respectively 

normal midplane strains in x and y directions . respectively 

6 rotation angle of the fibers with respect to the plate axis 

9’ rotation angle of the fibers in the inner laminas with respect 

to plate axis (See figures lU and 15.) « 

'^x**^y*'^xy plate curvature as defined by equation (2.5) 

A cross-flow angle 

X flutter parameter defined by equation (3.16) 

'’is poisson's ratio for orthotropic lamina defined by - ^2^®! 

®x* ®y* ^z stress components in x, y, and z directions, respectively 

^xy* ^xz* stress components 

u natural frequency 

ftindamental frequency of simply supported beam as given 
by equation ( 3 . 17 ) 


IC ,tc ,x 
X* y* xy 



LIST OF TABLES 


TABLE PAGE 

I. LAMINA MATERIAL PROPERTIES USED IN ANALYSIS T8 

II. MATERIAL PROPERTIES FOR COMPOSITE STIFFENED ALUMINUM PLATE. . 79 

III. MATERIAL PROPERTIES OF REFERENCE PLATE 80 

IV. COMPARISON OF FOUR LOWEST NATURAL FREQUENCIES FOR A SQUARE 

ANGLE-PLY PLATE °12^®22 ' *^12 ' 

V. COMPARISON OF FOUR LOWEST NATURAL FREQUENCIES FOR A 

SQUARE CROSS-PLY PLATE (Gi2^®22 * -5. ' *25) 82 


vll 



LIST OF FIGURES 


1 

I 


FIGURE f*AGE 

1. - Geometry of angle-ply plateu 83 

2. - Geometry of cross-ply plates 84 

3. - Lamina coordinate system 8^ 

4. - Lamina notation and positive direction of stress 86 

5. - Rotation of lamina axes 87 

6. - Positive direction of stress and moment resixltants 88 

7. - Plate geometry and coordinate system 89 

8. - Nomenclature and geometry of symmetric plates 90 

9. - Nomenclature and geometry of angle-ply plates 91 

10. - Nomenclature and geometry of general laminated plates .... 92 

11. - Nomenclature and geometry of composite stiffened aluminum 

plate 93 

12. - Convergence of solution for a symmetric, boron-epoxy, square 

plate with K = 4 95 

13. - Coalescence of frequencies for a symmetric, boron-epoxy, 

square plate with K = 4 96 

14. - Comparison of flutter solutions for square reference plate. 

98 


15. - Converged fluttei boundaries for square reference plate ... 99 

16. - Flutter boundaries for square, symmetric plates. N * N » 

N 100 

xy 

17. - Flutter boundaries for boron-epoxy, symmetric plate with 

a/b » 2.0. N « N if _ = 0 102 

X y xy 

18. - Flutter boundai;les for square, bcron-epoxy, symmetric plate 

with inplane normal loads. N = N »0;K*4 103 

y 


viil 



t 


lx 


C 

a 


FIGURE PAGE 

19 . - Flutter bound€U'ies for square, boron-epoxy, symmetric plate 

with inplane sheeu: loads, n »N ■0;K*4 lOU 

y X ’ 

20. - CoaJ-escence of frequencies for a symmetric, boron-epojy, 

square plate with inplane Jhear loads. N =N b 0;K»4; 
and 0 * 15 ° * 105 

21. - Flutter bovindaries for square, boron-epoxy, symmetric plate 

with crosn-flow. if *if =if »0;K*4 IO 8 

X y xy * 

22. - Flutter boundaries for square, angle-ply plate. 

if =: if = if =0 109 

X y xy 

23 . - Flutter boundaries for boron-epoxy, ongle-ply plate with 

a/b = 2.0. if = N * N * 0 Ill 

X y xy 

24. - Flutter boundaries for square, boron-epoxy, angle-ply plate 

with inplane loads. K = 4 112 

25 . - Flutter boundaries for square, boron-epoxy, angle-ply plate 

with cross-flow. N *N ®N =0;K = 4 Il4 

X y xy ’ 

26 . - Flutter bo;jndaries fc" squ6u*e, boron-epoxy, general 

laminated plates, ifsifeif =0 115 

X y xy 

27 . - Comparison of flutter boundaries for square, symmetric, 

angle-ply, and general laminated plates. N = N ■ N * 0 . 117 

X y xy 

28 . - Effects of cross-flow on flutter of square, symmetric, 

angle-ply and general laminated plates. N ■ N « N 0, , II 8 

X y 3QT 

29 . - Flutter boundaries for square, condos ite-stiffened aluminum 

plates. N = N * N ® 0; mass . . /mass , - 

X y xy aluminum con5>osite = 1.0 119 


Chapter 1 


IMTRODUCTION 


Laminated compoaite materials have been in use for many centuries. 
Pieces of laminated wood have been foimd which date back to about 1300 B.C. 
Some excellent swords were made in the 15th centviry by laminating several 
layers of steel to provide an extremely hard and keen cutting edge with 
a softer tough body. Although some of the advantages of using composite 
materials have been known for some time, only in recent years have deter- 
mined efforts been made to fully develop composite materials. The advent 
of high performance aircraft and spacecraft has brought about the need 
for more efficient use of materials. Recent energy shortages have em- 
phasized in even a greater way the Importance of developing composite 
materials for everyday usage. 

In recent years, composite materials have been the subject of a large 
number of experimental and analytical Investigations and have been develop- 
ed to the point where structural components are now being used in actual 
applications. Previously, most of the composite applications have been 
structural components built to replace existing conventional components. 
However, new designs are emerging which have considered composite materials 
in the total derign. As new applications ai'e considered for composite 
designs, a better understanding of the behavior of structural components 
under various load conditions is needed. 

One structm'al component that is used extensively in a wide variety 
of applications is the flat plate. Plates fabricated from composite 
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materials usually consist of individusd lamina bonded together. Each 
lamina usually consists of a number of high strength filaments adigned in 
the same direction and held in place by a plastic matrix materisd. 

The lamina may be highly orthotropic, and since the principsd material direc- 
tion or the angle the fibers make with the plate axis may differ for each 
lamina, the resulting anisotropic plate is more complex to analyze than 
an orthotropic plate. 

This study is being conducted to obtain a better understanding of 
the flutter characteristics of laminated composite flat plates. However, 
the following discussion will include a general review of previous studies 
of laminated composite plates. Most previous studies of laminated plates 
have been concerned only witn deflection, vibration, and buckling char- 
acteristics. The development and solution of the governing equations 
for plates under various conditions will be reviewed and a proposed area 
of investigation identified. Solution techniques and possible problem 

floreas will also be discussed. 

O 

0 1.1 Review of Literature 

Numerous theoretical analyses of composite laminated plates have been 
made and are discussed in references 1-23. These references represent 
work that cover the spectrum from classical snail deflection theory and 
approximate solutions to non-linear theory and exact solutions. Much of 
the early analytical work for composite plates was conducted using classi- 
cal small deflection theory, and solutions were obtained only for special 
classes of plates which resulted in considerable simplification of the 



analysis. Racent studies have obtained solutions to more general plate 
problems with transverse shear and rotary inertia effects. 

In addition to the reports and articles, some excellent books have 
considered laminated composite plates. Diet (2U) in 19^9 published a good 
review of much of the early analytical work in laminated composites and 
gave a thorough listing of references. Ashton, Halpin, and Petit (25) 
in 1969 published an excellent primer for those new in the field. They 
give a good development and discussion of the governing equations. Ashton 
and Whitney (26) in 1970 published a good review of the progr?ss in lami- 
nated plate theory and the various solutions and solution techniques. 

1.1.1 Formulation of Equations 

Two early investigators. Re j saner and Stavsky (l), analyzed a plate 
composed of two orthotropic laminas of equal thickness allg..ed so that the 
principal material axes of the two laminas are rotated at an angle of 
and -6 with the plate axes (later referred to as angle-ply plates). 

They found that even for a small deflection theory analysis, there exists 
a coupling between the transverse bending and inplane stretching equations 
that does not exist for orthotropic plates. They shoved that the coupling 
effect exists in the boundary conditions as well as in the governing 
equations. 

Relssner and Stavsky (l) formulated the governing equations as two 
Uth order partial differential equations in terms of the Airy stress 
function and the transverse displacement. Due to the coupling of the 
two equations, they must be solved simultaneously, and inplane, as well as 
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transverse Ijoundary conditions, must be specified at each boundary to 
obtain a particular plate solution. Stavsky (2) forrulated the govern- 
ing equations as an eighth order partial differentied equation in terms of 
the Airy stress function and later (3) in the form of three simultaneous 
partial differential equations in the displacements n, v, and w. Al- 
though each of the formulations have advantages in certain problems, 
they present about the same order of difficulty in obtaining a solution. 

Later contributions to the laminated plate formulations have been 
made by several authors. Tsai euid Azzi (U) added ' .ermal stresses to 
the formulation of laminated plate equations. Whitney and Leissa (5), 
and Yang, Norris, and Stavsky (6) included rotary inertia terms in the 
formulation amd the latter authors etlso included transverse shear which 
becomes important for thick plates. The formulation by Yang, Norris, and 
Stavsky results in five coupled partial differential equations in terms of 
the displacements u, v, and v and two rotations of the normals to the 
midplane. 

1.1.2 Solutions for Symmetric Plates 

Although available formulations of the governing equations and 
boundai*y conditions are applicable for general laminated plates, most 
of the solutions have been obtained for special classes of plates which 
result in considerable simplification to the governing equations. 

Stavsky (2) in his early work showed that for laminated plates symmetric 
about the geometric midplane both in properties emd fiber orientation, 
the transverse bending and inplane stretching equations and the boundary 



conditions become uncoupled and ccm be solved Independently. However, 
he shoved that even for symmetric laminated plates, the governing equations 
and boundeury conditions still contain some oss-stiffness terms due to 
the laminas which result in the governing equations being different 
from those of an orthotropic analysis. He further showed that only for 
veiy special kinds of symmetric plates do the cross-stiffness terms vanish 
and the governing equations become identical with the orthotropic equa- 
tions. As a result of these simplifications, symmetric plates have been 
studied extensively in the literature. 

The Raylelgh-Ritz method of solution with beam mode shape functions 
has been used effectively by several investigators to obtain deflection, 
vibration, and buckling solutions. Ashton and Waddoups (7) investigated 
symmetric laminated plates and presented analytical results for uniaxially 
loaded plates with the loaded edges clamped and the unloaded edges free. 
They also presented calculated freguencies and mode shapes for canti- 
levered plates which showed good agreement with experiment. Ashton (8) 
extended the emalysls to include nonuniform cross section and material 
properties and presented resvilts for tapered plates under imaxial 
lo ^ing with simple supported and clamped boundary conditions. Ashton 
(9) further extended the analysis to include elastically restrained 
boundary conditions and presented some deflection results with re- 
stredned boundary conditions. 

Later Srinivas wd Rao (lO) obtained closed form solutions for the 
free vibration and buckling of simply supported symmetric laminated 



thick plates using lineeu: small deflection theory. They compared their 
results with thin plate theory and showed that the error associated 
with thin plate theory Increases with plate thickness and the modular 
ratio between the lamina. A similar analysis by Whitney (11) showed 
the same trends. Thus the definition of "thin” for laminated plates 
must take into account modular ratio as well as plate thickness. 

1.1.3 Solutions for Angle-ply and Cross-ply Plates 

Angle-ply plates are unsymmetric and consist of an even nxanber of 
layers having the same thickness and elastic properties, and the ortho- 
tropic axes of symmetry for each ply are alternately oriented at angles 
of +0 and -6 to the plate axes (See fig. 1.). Cross-ply plates are un- 
symmetric and consist of an even number of layers all of the same thick- 
ness and elastic properties with the orthotropic axes of symmeti in 
each ply alternately oriented at 0 degrees and 90 degrees to the plate 
axes (See fig. 2j. Whitney and Leissa ( 5 ) showed that for angle-ply 
plates and cross-ply plates, the governing et^uations and boundary condi- 
tions are considerably simplified but SL.’e still coupled and must be solved 
simultaneously. As a result of the simplifications, euigle-ply and cross- 
ply laminated plates have been studied by several investigators to deter- 
mine the effects of bending-extensional coupling. 

Several solution procedures have been developed which are effective 
in solving the coupled governing equations and boundary conditions. 

Whitney and Leissa ( 5 ) and Whitney (12) used a Fourier series technique 
to obtain solutions for deflection, free vibration, and buckling of both 
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angle-ply and cross-ply plates for simply supported boundary conditions 
which allow inplane displacements normal and tangent to the boundary for 
cross-ply and angle-ply plates respectively. Bert and Mayberry (13) used 
a Rayleigh-Ritz approach with beam mode shape functions to obtain free 
vibration results for cross-ply and angle-ply plates with clamped 
boundary conditions. Whitney (lU) used a Galerkin procedure to obtain 
solutions for the shear buckling of cross-ply plates with simply suppor- 
ted boundary conditions. 

In all cases cited, the bending-extensional coupling significantly 
affected the results. Bending-extensional coupling has the overall 
effect of reducing the plate stiffness and thus increasing the static 
deflection and reducing the natural frequencies and buckling loads. The 
effect of coupling is increased as the number of laminas are decreased 
and as the degree of anisotropy between lamina is increased. 

Fortier (15) used a Rayleigh-Ritz solution procedure to investigate 
the effects that various types of inplare boundary restraints have on 
the behavior of angle-ply and cross-ply plates with small initial curva- 
ture. He also considered the effects of treuisverse shear and inplane 
loads on the natural frequencies. He found that inplane boundary condi- 
tions, initial curvature, transverse shear, and inplane loads all have 
a significant effect on the behavior of unsymmetric plates. However, 
Whitney (16) considered transverse shear for cylindricsd bending of 
symmetric and cross-ply pistes and showed that transverse shear has less 
effect on the deflection of cross-ply plates than symmetric plates. 
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l..j U Approximate Solutions 

Due to the difficulty of obtaining general plate solutions, approxi- 
mate solutions have emerged also. A "reduced bending stiffness" method 
was proposed by Ashton (20) in which the bending-extensional coupling 
is accounted for in an approximate way. Approximate solutions were ob- 
tained by reducing the bending stiffness of the plate by an amount deter- 
mined by the coupling terms and then neglecting the coupling effect. 
Ashton used the method to compare maximum deflections of a cross-ply and 
an angle-ply plate with measured results and obtained good agreement. 
Whitney (21) also used the method to compare results with those calcu- 
lated using the Fourier series technique and showed good agreement for 
5oth angle-ply and cross-ply plates for simply supported boundary condi- 
tions, but for lamped boundary conditions, relatively poor agreement was 

obtained for certain angle-ply cases. 

A more refined amedysis was presented by Srinivas (22) which con- 
siders transverse shear deformations and rotary inertia effects. By 
assuming piece- ds« linear variations of the displacements u and v 
and constant •''^alues of w across the plate thickness, the problem be- 
comes two dimensional. Trigonometric series solutions were obtained for 
the approximate two dimensional problem for plates with simply supported 
oundary conditions, and the results showed good agreement with exact 


theory . 
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1.1.^ Dynamic Stability 

A limi ed number of investigators have considered the dynamic 
stability of laminated plates. Smirnov (17) considered the flutter of 
an infinite sandwich plate subjected to cylindrical bending in a gas 
s\.ream and obtained an exact solution. He later extended the analysis 
(18) to include semi-infinite sandwich plates with veur*ious boundary con- 
ditions. He used linear piston theory aerodynamics and obtained solu- 
tions for both clamped and simply supported semi-infinite plates. 

Librescu and Badoiu (19) and Ramkumar (23) analyzed the flutter of 
flat, symmetrically laminated, simply supported plates with inplane 
normal loads. They used linear piston theory dynamics with aerodynamic 
damping and the Rayleigh-Ritz solution procedure with beam mode shape 
functions. Flutter boundaries are presented as a function of the 
orientation of the principal axis of the lamina for various inplane 
loads and aerodynamic damping coefficients. However, it should be 
noted that all the flutter boundaries presented in reference 19 were 
obtedned using only two terms in both the streainrise and cross-stream 
direction and the results presented in reference 23 were obtained using 
ten terms in the strearawise direction but only two terms in the cross- 
stream direction. Thus, in both cases, the results presented may not 
have been completely converged. 

1.2 Statement of Problem 

Althou^ much work has been done in the analysis of composite 
laminated plates, as discussed in the literature survey, panel flutter 



stands out as having received very little attention. Panel flutter has 
long been recognized as a problem for the design of conventional, panels 
subjected to supersonic flow. This is evident by the laurge number of 
reported panel flutter investigations some of which are listed in refer- 
ences 27-U2. The panel flutter work that has been done for laminated 
plates is so limited that it is of little value to the designer. It is 
evident that to design efficient laminated con^osite panels for super- 
sonic application, a better understanding of their flxitter characteris- 
tics is needed. 

Another area that has received little attention is the solution of 
general laminated plate problems wher^ the plates may have any nimiber 
of layers stacked in an surbitrary sequence with the fibers in each layer 
rotated at an arbitrary angle to the plate axis. For general laminated 
plates, the complete bending and extensional governing equations and 
boundary conditions are coupled and must be solved simultaneously to 
obtain plate solutions. Considerable difficulty is encoiintired in 
solving the coupled equations emd boiuidary conditions. Thus, edthou^ 
the governing equations and boundary conditions have been formulated for 
generetl laminated plates, no general solutions were found in the litera- 
tvire survey. All the solutions presented were for specially laminated 
(symmetric, angle-ply, or cross-ply) plates whose special geometry 
results in simplifications to the governing equations and boundary con- 
ditions which make it considerably easier to obtcdn solutions. However, 
since in the design of laminated plates it is not always practical or 
desirable to use specially laminated plates, an suialysis is needed which 
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obtains solutions for general laminated plates. 

The purpose of this investigation is to obtain the flutter charac- 
teristics of flat, general laminated plates using approximate methods. 

The analysis will be based upon small deflection theory but will consider 
the general coupled governing equations and boundary conditions for 
simply supported plates. Thiis, an approximate solution procedure will 
be developed for completely general plates which may have any number of 
layers stacked in an surbitrary sequence and with the fibers in each 
lamina rotated at an arbitrary angle to the plate sixes. The p]ate may 
be subjected to uniform inplane normal and shear loads. 

An extended Galerkin method will be used to obtain approximate 
solutions to the governing equations and boundary conditions. The 
extended Galerkin method will be used since it provides a straight- 
forward solution procedure for nonconservative problems (See ref. U 2 O 
using simple series to describe the displacements. The aerodynamic 
loading on the panel will be assumed to be given by linear piston theory 
and the flow may be at an arbitrary cross-flow angle. Piston theory 
aerodynamics will be used because it gives simple expressions for the 
aerodynamic loading and has been shown to be applicable for Mach numbers 
greater thaui 1.6 (See ref. 28 and 29») . 

Since symmetric and angle-ply plates have been the subject of many 
investigations, approximate solutions will be obtained for typical 
laminated plates from each class, and a limited parametric study will be 
conducted to determine the effects on the flutter boundaries. The param- 
eters to be studied will include the number and orientation of the plies, 


length-width ratio, inplane shear and normal loads, and cross-flow 
angles. The flutter boundaries from this analysis will be con^pared with 
those calculated using classical orthotropic plate theory and the reduced 
bending stiffness method. 

Since general laminated plates have not been investigated prior to 
this time, approximate solutions of their governing equations will be of 
special interest. Thus, flutter boundaries will be calculated for some 
typical general laminated plates and the results will be congjared with 
those obtained for specially constructed symmetric and angle-ply plates. 
Flutter boundaries will also be calculated for some composite stiffened 
aliuninxim plates that do not have a syijnetric or angle-ply construction. 

Althou^ the approximate procedure will be developed for the purpose 
of obtaining flutter solutions, the analysis will also have the capability 
of giving natural vibration frequencies, and inplane normal and shear 
static buckling loads. Since no numerical resxilts are available in the 
literature for general laminated plates, a significant contribution to 
the literature could be made by using the natural vibration and static 
buckling capabilities of the analysis. However, only a limited number 
of natural frequency calcvilations will be made to compare with published 
results to verify the solution procedure. 
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Chapter 2 

DEVELOPMENT OF GOVERNING EQUATIONS 

In order to develop the governing equations for this investigation, 
several assumptions are msuie as follows: 

1. The plate is constructed of flat, uniform thichness layers of 
orthotropic sheets bonded together. The direction of principal stiffness 
of the individual layers do not in general coincide with the plate edges. 

2. The plate is thin; i.e., the thickness is much smaller than the 
other physical dimensions. 

3. The displacements are small compared to the thickness. 

U. Each lamina obeys Hooke's law. 

5. The Kirchhoff hypothesis is used; i.e., normals to the midplane 
of the undeformed plate remain straight and normal to the midplane dxurlng 
deformation. 

6. Transverse shear and normal strains are negligible. 

7. Body and rotary inertia forces are negligible. 

8. The plate is of constant thickness. 

These assumptions give rise to the conclusions that y i » Y , x , and 

xz yz xz yz 

a are negligible which is the case for an approxlsiate state of plane 
stress. 

The coordinate system used to Identify the plate and lamina 
geometry is shown in figure 3. The distances to the Individual lamlnas 
are measured from the geometrical midplane of the plate. The positive 
directions for stresses are shown in figure U. 
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2.1 Lamina Stres8>Strain Equations 


From elementary strength of materials, Hooke's law for an orthotro- 
pic lamina in a state of plane stress is given as follows (See ref. 
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where the stresses ®i*®2*^12* strains ^12 

referred to the direction of principal stiffness and the subscript J 
refers to a particular lamina. This relationship, written for another 
(x-y-z) system of axes (See figure 5.) where the x-y axes are rotated at 

an angle 0 with respect to the 1-2 axes, is given as follows: 
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( 2 . 2 ) 


where the matrix is a transformed matrix with the elastic proper- 

ties along the principal axes of the lamina rotated to the x-y-z 
(plate) axes system. The matrix has non-zero values for all terms 
when the principal axes of the lamina do not coincide with the plate 
axes system (for example, a fibrous composite with the fibers rotated at 
an angle with the plate axes). 
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From elementary strength of materials considerations, the strain at 
any point in a lamina undergoing deformation can he expressed in terns 
of the deformation of the geometrical midplane of the lamina. For 
small deflections and a state of plane stress, the strain at any point 
in a lamina z distance from the midplane is given in terms of the 
midplane strain and curvature as follows: 
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where e e and y ® ere the midplane strains and < »< » and k 

DC 3C 

are the midplane curvatures. The midplane strains and curvatures de- 
fined in terms of the displacements are needed in the later development 
of the equations and are given as follows: 
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Using equations (2.2) and (2.3)* the stress state at any point in a 
Isnlna SMy be written in terms of the midplane strain and curvature as: 
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( 2 . 6 ) 


or 






(2.7) 
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where z is the distance from the midplane to the point. 

2.2 Lamina and Plate Constitutive Equations 

Since for a plate cosqposed of several lemlnas, the stresses, 
strains, and displacanents will be different for each lamina, it is 
convenient to define some equivalent system of forces and moments which 
will be considered to be applied to the midplane of the plate. Within 
the plane stress assumption, such a system will he defined for a lamina 
in terms of three stress and three moment resultants. The stress and moment 
resultants for a lamina are defined as: 

(«) 




( 2 . 8 ) 
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(b) 


(c) 


(•) 


(b) 
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(2.9) 


and arc statically equivalent to the actual stresses on the lamina and 
may be considered to be applied at the midplane of the lamina. The 
positive directions for the stress and moment resultants on the plate 
are shown in figure 6. Using equations (2.6), (2.8) and (2.9) the 
stress and moment resultants can be related to the strain and curvature 
of the midplane. The resulting relations are known as the lamina 
constitutive equations and are given as follows: 
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The plate constitutive equations will be developed from equations (2.10) 
and I (2.11) by sunnaing up the individual terms for each lamina. 

Thus» 
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since e*’ and K are functions of x and y only and Q is a function of 
lamina properties, they can he taken outside the integred. Thus, 
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or after Integrating these may he vritten aa: 
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Writing the constitutive equations for laminated plates in matrix 
notation, they become: 

[N] a [A] [e®] ♦ [B] [<] (2.19) 

[M] a [B] [e®] + [D] [ic] (2.20) 

These equations indicate that for a general laminated plate either stress 
or moment resultants produce both strains and curvatures. Thus there is 
a coupling between the stress end moment resultants through the matrix [B]. 
Since [b] Is an even function of h (Saa aq. (2.l8),), it is zaro and hence 
the equations are uncoupled for classical orthotropic plates and for 
laminated plates that are symmetric about the plate midplane. 

Equations (2.19) and (2.20) will be rewritten in a more useful form 
for the following development by solving equation (2.19) Tor the midplane 

strains, and rewriting equation (2.20) in terms of the stress resulteuits. 
Equations (2.19) and (2.20) become: 
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[£°] - [A*] [H] + [B*] [ic] (2.21) 

[M] = [N] + [D*] [ k ] (2.22) 

where 

[A*] = [A“^] 

[B*l = -[a'^1 [b] 

[D»] = [D] - [B] [a'^1 [B] (2.23) 

2.3 Plate Governing Equations 

Due to the coupling between moment euid stress resulteuits noted in 
equations (2.19) and (2.20), two differential equations are required to 
describe the behavior of the plate \uder various load conditions. The 
first equation will be developed using plate equations of motion and the 
second equation will be developed from the compatibility conditions that 
must be satisflea between the strains. 

2.3.1 Equations of Motion 

The equations of motion for a thin laminated plate are identical to 
those of homogenous plate theory. For a plate -subjected to inplane normal 
and shear loads and loads perpendiculeu: to the plate surface, the .equations 
of motion are given as follows (See ref. 26^ where the inertia terms in the 
X and y directions are considered to be negligible. 
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The inplane loads N . N • and N may he separated into the uni- 

X y xy 

form externally applied loads arid the loads induced by the plate deflections 
as follows; 
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Substituting the expressions for N , N , and N given by equation (2.2?) 

X y 3QT 

into equations (2.2U), (2.25), and (2.26) results in the following 
equations of motion: 
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The higher order terms 
have been neglected in equation (2.30). 

Defining a stress function F such that: 
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then equations (2.28) and (2.29) are satisfied identically by the function. 
The function F is referred to in the literature (See ref. 26.) as the 
Airy stress function. Using tbe definitions of F given by equations 
(2.31), the definitions of cuivature given by equation (2.5), and the 
expression for the moments gjiven by equation (2.22) in equation (2.30), 
the following expression is obtained: 
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If p is taken as the aerodynamic loading due to air flow over the 
plate surface* then equation (2.33) is one of the equations needed to 
solve for the flutter of flat laminated plates. Aerodynamic loads pre- 
dicted by piston theory have often been used in the literature (See ref. 
28 - l+2j to obtain flutter solutions. Piston theory gives a relatively 
simple expression for p and has been shown to be applicable for Mach 

numbers greater than 1.6 (See ref. 28 and 29.). For supersonic flow 
over one side of a flat plate at an emgle of A with the x-axis as 



shown in figure 7» piston theory gives the following expression for the 
aerodynamic loads. 
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Substituting the expression for p given by equation {2.3h) into 
equation (2.33), the loilowing governing equation is obtained for flutter 
of flat laminated plates. 
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This expression is a l*th order partieil differential equation which governs 
the flutter behavior of flat laminated plates but is in terms of two un- 
knowns w and F. Thus em additional equation relating w and F is 
needed before solutions to flutter problems may be obtained. 
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2.3 >2 Compatibility Equation 

The additional equation needed to solve for the plate flutter 
behavior is obtained from the compatibility condition that must be 
satisfied between the midplane strains. The appropriate equation is 
obtained by eliminating the midplane displacements from the strain- 
displacement relations given in equation (2.4). Differentiating 
(eq. 2.4(a)) twice with respect to y and c° (eq. 2.4(b)) twice with 

y 

respect to x and adding, the following expression is obtained: 
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Equation (2.37) is referred to as the compatibility equation and will be 
used to develop the second governing equation in terms of w and E. 
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Using the expressions for strains given hy equation (2.21) and the 
force resultants given hy equations (2.27) and (2.31) and the 
ctirvatures .■ given hy equation (2.5) In equation (2.37) • give rise to the 
following expression: 
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This Is the second governing equation needed to solve for the plate be- 
havior. It Is also a fourth order partlsQ. differential equation In terms 
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of F and w and for flutter problems must be solved simultaneously with 
equation (2.35). In order to solve the equations for a specific problem 
the boundeuy conditions must be specified. 

2.U Boundary Conditions 

Proper boundary conditions which guarantee \mique solutions to the 
two governing equations must be specified. It has been found (See 
reference 26.) that the necessary boundary conditions are those of 
classical homogeneous plate theory plus those of an inplane elasticity 
problem. Thus due to the coupling, the usual deflection, slope, moment, or 
shear boundary conditions used in classifical plate theory do not give 
unique solutions but inplane boundary conditions must also be specified. 

The resulting boundary conditions require one member of each pair of the 
following quauititles to be specified along the boundaries: 
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where n and t are used to denote coordinates normal and tangential to the 


plate boundary respectively. 
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For this analysis, the boundary conditions chosen represent those 
for a simply supported plate with no inplane edge restraints and are 
given as follows: 


V (O.y) • w (a,y) ■ >.(x,0) • w(x,b) ■ 0 (a) 

(O.y) - (a.y) • 0 (b) 

M (x.O) • M (x.b) • 0 (c) 

(O.y) « (a.y) • 0 (d) 

(x.O) » r (x.b) . 0 (e) 


N' (o.y) » N' (a.y) * N' (x.O) ■ N’ (x.b) ■ 0 (f) 

AJ 3QT AgT 


(2.1»1) 


These equations will be rewritten in terms of w and P using the expression 
for momentr given by equation (2.22) and the definitions of strain, cur- 
vature. and F given by eqiuations (2. it), (2.3). and (2.31) respectively. 
Thus the boundary conditions become: 

w (o.y) s w (a.y) ■ w (x.O) « w(x,b) • 0 (2.U2) 
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The governing equations (2.35) «nd (2,39) along with the boundary con- 
ditions (2.^2) through (2.U7) will be solved by using approximate techniques 
in the next chapter. 
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Chapter 3 

APPROXIMAT F SOLUTIONS OF THE FLUTTER EQUATIONS 

The two governing fourth order partial differential equations (2.35) 
and (2.39) along with boundary conditions given by equations (2.U2) 
through (2.1*7) will be used to obtain approximate flutter solutions. The 
analysis will be developed for the flutter of general laminated plates where 
erbitrary stacking sequence and orientation of the lamina fibers is 
permissible. Since certain special classes of plates have been studied 
almost exclusively in the literature, the simplifications to the 
general solution due to the special constructions will be discusse-. The 
special classes of plates considered are symmetric laminated plates and 
angle-ply plates. Symmetric plates may be composed of any even munber of 
lamina of arbitrary thickness and orientation of the fiber directions as 
long as for every lamina above the plate midplane there is an identical 
(in thickness, materj.al properties, and orientation) lamina equal 
distance below the plate midplane. Angle-ply plates are less general than 
symmetric plates and consist of an even number of laminas all of the same 
thickness and elastic properties with the orthotropic axes of eyinBietry 
in each lamina alternately oriented at -*-6 and -6 to the plate axis. 

Although the analysis and solution procedure will be developed basically 
to study panel flutter, the anadysis is general ar.d the result^'ng computer 
program can also be used to determine static buckling loads and natural 
vibration frequencies. 
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3.1 General Laminated Plates 

An extended Galerkin u»thod will be used to obtain approximate 
flutter solutions since it provides a strai^t-forward solution proce- 
dure for nonconservative problems (See ref. U2.) using a simple series 
to describe the assumed displacements. The extended Galerkin method 
admits solutions only when the assumed deflections satisfy the geometric 
(deflection and slope) boxondary conditions. However, the ntimber of 
terms required in the solution to obtain converged results are probably 
reduced if the assumed deflections satisfy some or all of the natural 
(force or natural constraints) boundary conditions. Since functions are 
not available which satisfy all the boundary conditions given by equa- 
tions (2.42) through (2.17), the capability to accoiint for the natural 
boundary conditions makes the extended Galerkin method particularly 
suited to this analysis. 

extended Galerkin method of solution is illxistrated by the 
following equation written in terms of the virtual work for displace- 
ments w. 



(governing equation )fiv dxdy ♦ 
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The terms under the double and single Integrals are usually Individually 
set equsd to zero and referred to as the governing equation and boundary 
conditions, respectively. For the boundary conditions to be satisfied, 
either the functions f.(w,x,y) or the variation of the deflections or the 
slopes must be zero at the appropriate boundaries. In usual applications 
of the Galerkin method, a series of functions is chosen which satisfy all 
the boundary conditions and the coefficients of the functions are deter- 
mined so that the double integral term is zero or equivalently the coeffi- 
cients are determined so that each term of the series is orthognal to the 
exact solution (See ref. 1»3.). When all the natural boundary conditions 
are not satisfied, the extended Galerkin method requires that the coeffi- 
cients be determined using the unsatisfied single integral terms as well 
as the double integred terms in equation 3.1* 

For the flutter of general laminated plates, a series of functions 
for both w and F are needed which exactly satisfy the geometrical 
boundary conditions given by equation (2.U2) and which satisfy as closely 
as possible the natural boundary conditions given by equations (2.U3) 
through (2.U7). The functions assumed for this analysis are given as 
follows : 
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These functions exactly satisfy the geometrical houndeury conditions and 
the natural boundary conditions given by equation (2.147) but the re- 
maining natural boundary conditions are not completely satisfied. The 
unsatisfied natural boundary condition terms are given as follows: 
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The boundary condition terms that must be included in the extended 
Oalerkin method written in terms of the virtual work are given as follows; 
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where the displacements u® end v® are given as follows and were obtained by 
substituting the definition of curvature given by equation (2.5) and the 
expression for strain given by equation (2.22) into equation (2.U) and 
integrating. 
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Using the assumed functions given by equations (3.2) and (3.3) and applying 
the extended Galerkin method, the orthogonality relations for the governing 
equations (3.35) (3.39) with the boundary condition terms given by 

equations (3.5) through (3.10) added may be written as follows; 
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Equations (3.1^) and (3.15) results in 2(MxN) equations in terms of the 

two unknowns C and H . In order to get the equations in a form 
rs rs o ^ 

that can be solved, H must be solved for in terms of C . This can 

rs rs 

be done using equation (3.15) and writing it in the following matrix 
form: 


[b ] [H ] = [a ] [C ] 
mnrs rs mnrs '■ rs 


(3.18) 


where r and s are summed from 1 to M and N respectively and 

b and a eire the coefficient matrices for H and C res- 

mnrs nmrs mn mn 

pectively. To obtain as a function of C^g» premultiply each side 

of equation (3.18) by the inverse of [b ] as follows: 
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Equation (3.20) results in a series solution for H in terms of the 

mn 


unknown coefficients C . 
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Using equation (3.20) to eliminate H from equation (3.1^), MxN 

mn 

homogenous equations are obtained in terms of the unknown The 

resulting equations may be written in the following form: 



where the c are the coefficients of C . Equation (3.21) is recog- 
mn mn 

nized as the well known characteristic equation emd has a non-trivlal 
solution obtained by setting the determinant of the coefficient matrix 
to zero. This results in an eO-gebraic equation of order MxN in terms 
of w. The roots of the algebraic equation sure the characteristic values 
or eigenvalves and represent the plate nondimens ionsd vibration fre- 
quencies . 

Flutter solutions have been obtained using a digital computer to 
solve equation (3.21) for its eigenvalves. Since the coefficients, 

Cj^j , sure a f\inctio:i of the flutter psurameter, X, the pla+e vibration 
frequencies are also a function of X. The criteria used to define the 
point of flutter is the lowest value of X which results in two of the 
frequencies coalescing and thus having a negative imaginary 

value. This criteria. was selected because a negative imaginary fre- 
quency results in the assumed deflection function (See equation (3.2).) 
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4 I vt* 

becoming unbounded since it is multiplied by e^ . Similar criteria 
has been used extensively in the literature (See ref. 39 •) to define the 
point of flutter. 


3.2 Symmetric Laminated Plates 


Plates that are lami..ated such that they au*e symmetric about the 
midpleuie represent the largest class of laminated plates of practical 
use. Their symmetric construction results in warp- free structural 
elements which are desirable in most applications. Also, the governing 
differential equations and the boundary conditions are considerably 
simpler than for general lw.<ioated plates. Th\is, symmetric plates have 
been investigated mor^ any other class of laminated plates. 

As was pointed out in section (?.2) and as shown by equation 

(2.18 (b)), B is an even function of the lamina thickness and for 

mn 

symmetric plates all the terms vanish. Thus, the flutter equations 

(equations (2.35) and (2.39)) are considerably shortened and more impor- 
tantly become uncoupled. Then only equation (2.35) needs to be solved 
to obtain a flutter solution. For a symmetrically laminated plate, the 
governing flutter equation becomes; 
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where the asterisks have been left off the D terms since D « 

mn mn 

for symmetric laminated plates. 

Althou^ equation (3>22) is considerably sln^ler than the flutter 
equations for general laminated plates, it differs from the governing 
equation for flutter of an orthotropic plate by the a.lditional and 

Dgg terms. These terms appear in the governing equation as a result of 
the fibers not being aligned with the plate axes. 

Since for symmetric plates the governing eqxiations become uncoupled 
the boundary conditions that must be satisfied for a flutter solution 
are given by equations (2.U2) throu^ (2.hl«) and are rewritten as 
follows with the B terms removed: 
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These equations differ from the boundary conditions for a simply 
supported orthotropic plate by the ^16 eind ^26 terms. 

To obtain a stcmdard Galerkln solution for symmetric plates, it is 
necesscnry to assume a function for w that satisfies all the boundary 
conditions (3.23) through (3.25). However, the and D^g terms 

mal:e it difficult to .find sacisfactory functions. Therefore, the same 
deflection function used in section (3«l) for general laminated plates 
and given by equation (3.2) will be assumed Guid the extended Galerkln 
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method used to account for the unsatisfied natural boundary conditions. 
The boundary condition terms that must be accounted for by the extended 
Galerkin method are given as follows: 
b 


f o n ^ ^ A — 
j l6 9x3y 3x 


dy 


(3.26) 


a 

/ 


^ ^26 Ixly ^ ^ 


dx 


(3.27) 


Using the assiimed functions given by equation (3.2) and applying the 
extended Galerkin method to the governing equation (3.22) with the "16 
and boundary condition terms included, results in the following 

equation: 
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Equation (3>29) represents a system of (MxN) linear homogenous equations 
In terms of the unknown and Is recognized as the eigenvalue equa- 

tion. Flutter solutions are obtained for equation (3.29) in the same 
manner as discussed in section (3.1) for general laminated plates. 

3.3 Angle-ply Plates 

Angle-ply plates by definition are constructed so as to have an even 
number of layers ell of the same thickness and elastic properties an^ 
with the orchotropic axis of symmetry in each ply alternately oriented 
at +6 fiid -6 to the plate axis (See Fig. 2.). Although angle-ply 
plates do not represent as Isu-ge or as important a group as the symmetric 
plates, they are in?)ortant from an analytic viewpoint for obtaining a 
better understanding of general laminated plates. Because of the special 
geometry of ang] e-ply plates, the governing equations are considerably 
simplified but are still coupled. Thus, angle-ply plates retain many of 
the characteristics of general laminated plates but relatively simple 
solutions can be obtedned. These characteristics have resulted in the 
angle-ply plates being studied extensively in the literature to cbtain a 
better understanding of general laminated plates. 

Fur angle-ply plates, it can be shown (See ref. 5.) that this 
special construction results in some of the coefficients being zero as 
follows : 
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Setting these terms to zero in the flutter equations (2.35) and (2.39), 
resiilts in the following governing equations. 
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These e 4 uations are coupled and must he solved simultaneously. 

Since the governing equations for angle-ply plates are coupled, all 
the boundeu*y conditions given by equations (2.42) through (2.47' must be 
satisfied. However with some of the coefficients being zero as indicated 
by equations (3.29), the boundary conditions are simplified as follows: 
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If the functions for w and F used for the general laaiinated 
plates (section 3.1) and given by equations (3.2) euad (3.3) are used, 
then boundary conditions (3.33), (3.3^), (3.35), and (3.38) are satis- 
fied. Boundary conditions (3.36) and (3.37) are not satisfied directly 
but will be shown to be satisfied indirectly by the following develop- 
ment. Since equations (3.3b) and (3*37) are not satisfied, the terms 
that must be included in the extended Galerkin method are given by 
eqmtions (3.7) and (3.8) and are repeated as follows; 
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Using equations (3.9)» (3.10), (3.29) and (3.30), the displacements u° 
and v*^ are given as follows: 
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Using the assumed functions given by equations (3.2) and (3*3) in equa- 
tions (3.39) and (3.1+0) results in zero values for u° and v° on the 
boundaries x = 0 and a, and y = 0 and b, respectively. Thus, aJ.1 
the boundary conditions are satisfied without using the extended 
Galerkin method. 

Substituting the expressions for F and w given by equations 
(3.2) and (3.3) into equations (3.31) and (3.32) results in the follow- 
ing equations: 
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Eliminating between the two equations, the following expression is 

obtained: 
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Applying the Galerkin procedure to equation (3.^3) and reeirranging the 
terms, the Tollowing equation is obtained: 
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Equation (3.^5) results in MxN homogenous equations in terms of 
and is recognized as the eigenvalue equation. Flutter solutions are 
obtained for equation (3.^5) in the same manner as discussed in section 
3.1 for general and symaetric laminated plates. 


3A Approximate Solutions 


An approximate theory (See ref. 20.) referred to as the reduced 
bending stiffness theory will be compared with the present analysis. In 
the approximate theory, the plate bending stiffness is reduced by an 
amount depending on the coupling between the governing equations and 
then the coupling is neglected when solving the equations. The reduced 
bending stiffness is the same as that defined by equation (2.18(c)). 
Thus, the governing equation for the reduced bending stiffness theory 
is obtained from equation (2.35) by neglecting the coupling 
is given as follows: • 
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Since the coupling is neglected, only boundazy conditions given by equa- 
tions (2.42) through (2,44) with [B ] = 0 need to be satisfied. The 

inn 

bcundeuy conditions are given as follows: 

w(0,y) = w(a.y) = w(x,0) = w(x,b) = 0 (3.47) 
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The resulting governing equation (3.46) and boundary conditions 
(3.47) throu^ (3.49) are identical to those used for sysmetric lamin- 
ated plates (See equations (2.22) through (2.25).) except that [D*] Is 
used instead of [D]. Thus, the approximate theory reduces the solution 
for general laminated plates to that of a symmetric plate with the bend- 
ing stiffness reduced. For angle-ply plates; the “d ^*26 

(See equation (3.29).) are zero emd the approximate governing equation 
(3.46) and boundary conditions (3.47) through (3.49) are identical to 
those for classical plate theory. Thus, the reduced bending stiffness 
theory can be used to obtain solutions for angle-ply plates from 



piiblished classical plate theory resvilts. Since for symnetric plates 
[B^] - 0 and [D*] = [D], the reduced bending stiffness theory and the 
present analysis are identical. 

The advantages of using the approximate theory are evident for 
angle-ply plates since results can be obtained using well known published 
solutions. For symmetric and general laminated plates, the advantages 
are not so clear since a symmetric analysis similar to the one discussed 
in section 3.2 is necessary even to obtain approximate solutions. How- 
ever, prior to the present analysis the reduced bending stiffness theory 
was the only solution procedure available for general laminated plates. 
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RESULTS AND DISCUSSIONS 


The eqiiations governing the behavior of general laminated plates 
have been programed for approximate solution on a digiteJ. computer. 
Although the computer prograjns were developed largely to obtain flutter 
results, natural vibration frequencies and static buckling loads may 
readily be obtained using the programs* Special euialysis and programs 
were developed to obtain solutions for symmetric and angle-ply plates* 
However, the advantages of using the specieuL programs over the one 
developed for the geneial laminated plates were minimal and in all cases 
gave identiced results. 

Flutter boundaries have been calculated for symmetric, angle-ply, 
and general laminated plates. The plate nomenclature and geometry are 
shown in figure 8 for symmetric plates, in figure 9 for e^gle-ply plates, 
and in figure 10 for genered laminated plates. Flutter ooundeuries were 

also obtained for an aluminum plate with one, two and four layers of 

/* 

composite material applied to one or both sides of the plate as shewn in 
figure 3J-. The material properties of the individual lamina used in 
making the calculations cure typical of those for boron-epoxy and glass- 
epoxy materials. These material properties are given in table I. 

Material properties for the composite stiffened aluminum plate are given 
in table II. 


52 



U.l Convergence of Results 


When using a series solution, it is important to determine that 
sufficient terms are used in the analysis to obtain converged results. 
However, xising more terms than necessary results in considerably longer 
computation times. The usual procedure for determining convergence is 
to start with a small number of terms and increase the number until the 
solution does not change as more terms are included. For solutions 
obtained herein, a converged result is one which changes no more than 
1 percent when more terms are included in the flutter determinant. This 
method of determining convergence is demonstrated in figure 12 for a 
four-ply, symmetric, square panel where the flutter pareuneter is shown 
as a function of the number of terms in the x-direction. Curves are 
shown for 0=0° and 6o°, and the symbols represent the different 
numbers of terms used in the y-direction. For 0=0°, four terms in 
the x-direction and only one term in the y-direction are sufficient for 
converged results. For 0 = 6o°, approximately six terns in both the 
X- and y-directions are needed to obtain results converged within 1 per- 
cent. Note that for 0 » 60°, the converged solution is approached from 
higgler values of X for m > U and n > 2. This siiggests that the 
cross-stiffness terms have an adverse effect on the flutter solution, a 
phenomena that will be discussed in a later section. 

The modes that coalesce to produce flutter for 0=0° and 0 = 60° 
are shown in figures *13(a) and 13(b), respectively, where X is shown 
as a function of the non-dimensional frequency squared. For both cases. 



the (l»l) and (2,1) modes coalesce and result in flutter. The different 
nunibers of terms required to obtain converged results for the tvo cases 
are due to theD^g and D^g cross-stiffness terms. For 0 = 0°, the cross- 
stiffness terms are zero, and the plate Is an orthotropic plate for 
which It has been shown (See ref. Uo.) that only the first mode In the 
cross-stream direction has an effect on the flutter solutions. For 
0 = 60°, the cross-stiffness terms couple the modes In the x- and 
y-dlrectlons, and consequently, all modes effect the flutter solutions. 

An attempt has been made throughout the remainder of this study to in- 
sure that the values presented are converged within approximately 1 
percent . 

k,2 Con^>arlson with Literature 

In order to verify, as well as possible, the accuracy of the present 
analysis, the analysis will be \ised to obtain flutter data coiQ)arable 
with that presented in reference 19 for square symmetric plates. The 
material properties for the plate considered in reference 19 are given 
in table III. A comparison of the flutter boundaries is shown in figure 
lU where X is shown as a function of the orientation angle of the out- 
side lamina. The dashed curves were obtained from reference 19, and the 
solid curves were obtained using the present analysis. Both sets of 
curves were calculated using only two terms in the x- and y-directions. 
The agreement between the two analyses is very good except for 0* » 0° 
and 0* » 90° and for 0 near 90° where large discrepancies do exist. 
Since the present ansdysis uses a Galerkin method of solution and the 
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aruiysis of reference 19 uses a Raylei^-Ritz method, the discrepancy is 
attributed to the difference between the two methods when an insufficient 
number of terms are used in the euialysis to obtain converged results. 

For converged results, both analyses should give the same values. Con- 
verged flutter boundauries obtained using the present analysis are shown 
in figure 15. Although it is evident that the two-term solutions are 
not converged, there is still good eigreement between the trends shown by 
the two-term and the converged analysis. For 0=0° or 9 = 90° and 
for 0'= 0° or 0' = 90°, the present analysis should give the same re- 
sults as an orthotropic analysis since the cross-stiffness terns are 
zero. This is shown in figure 15 where the symbols represent orthotropic 
results obtained from reference 31. 

Although no flutter solutions were found in the literature for 
angle-ply or general laminated plates, some comparisons of natural fre- 
quencies for angle-ply and cross-ply plates are possible. The four 
lowest natural frequencies calculated using the present analysis are 
compared in tables IV and V with those calcxilated in reference 15 using 
a Rayleigh-Ritz method for square angle-ply and cross-ply plates, res- 
pectively. The natural frequencies calculated by the present analysis 
for simple supports should compare directly with those referred to in 
reference 15 as having hinge-free-tangential euid hinge-free-normal 
boundary conditions, respectively, for the angle-ply end cross-ply 
plates. For angle-ply plates, natural frequencies are given for the 
principal material directions rotated at 0°, 10°, 20°, 30°, and U0° 

with respect to the plate axes, llie agreement between the two analyses 
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Is excellent. 

For cro83-ply plates, natural frequencies calculated using the two 
analyses are compared (See table V.), for vedues of * 1, 10, 20, 

30, and ^0. The agreement between the two analyses is excellent except 
for the higher frequencies and hi^ values of the ratio E3^2^®22* 
ever, the differences sure too small to be of practical significance. 

Since a s'.ecial analysis was not made for cross-ply plates, the analysis 
for general laminated plates was \ised to obtain the natural frequencies. 
Thus, the assumed mode shape functions do not satisfy the boundary con- 
ditions exactly, but they are satisfied through the extended Galerkin 
procedure as discussed in section (3.1). The good agreement between the 
natural frequencies indicates the validity of lislng the Gederkin method 
to account for the boundary conditions not satisfied by the assumed mode 
functions. 

It. 3 Symmetric Plates 

Flutter boimdscries for synDetric plates with lamina properties 
typical of boron-epo^v materisil or a glass-epoj^ material sire presented 
in this section where X is shown sis a function of the angle that the 
fibers make with the x-cucis. The effects of inplane noiu.'. smd shear 
loads and cross-flow on the flutter of a squsure boron-epoj^ plate com- 
posed of ^ laminas are also shown. The values of X in each case sure 
based on ^rtiich is the value of with all the fibers aligned 


with the x-sucis 



1 


4.3.1 Flutter Boundaries 

Flutter boundaries for square, symmetric, glass-epoxy and boron- 
epoxy plates are presented in figure l6(a), and 16(b), respectively, for 
plates composed of two, four, and six laminae. For square symmetric 
plates, the hipest flutter stability is obtained with the fibers aligned 
with tho X-axis; rotating the fibers away from the x-axis results in a 
continuous reduction in flutter stability for values of 6 up to 90°. 
Flutter boundaries calculated using cleissical orthotropic plate theory 
which neglects the cross-stiffness terms and Dgg are seen to be 

Inaccxirate and nonconservative. Thus, inclusion of the cross-stiffness 
terms in the analysis has a destabilizing effect on flutter. Increasing 
the number of laminas composing the plate resvilts in improved flutter 
stability and flutter boundaries that are closer to those given by 
classlceQ. plate theory. P comparison of the flutter boundaries for the 
glass-epojQT (See fig. l6(a).> and boron-epoxy (See fig. 16(b).) plates 
shows that increasing the orthotropy ( for glass-epoxy material ^21^^22 

« 3, and for boron-epoxy material ®2i/^22 * °^ laminas causes 

the discrepancy between clMslcal plate theory and the present analysis 
to be more pronounced. Also, for the more hi^ly orthotropic material 
properties, orienting the fibers transverse to the flow has a larger 
destabilizing effect on flutter. 


Flutter boundaries similar to those presented in figure l6 for a 


square plate are presented in figxire 17 for a boron-epojgr plate with a 
length-width ratio of 2.0. These flutter boundaries show the solutions 
based i^on classical plate theory to be hi^ly nonconservative due to 
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the neglect of the crosi-stiffness terns. For the aA) “ 2.0 plate, 
rotation of the fibers avay from the x-axla resiilts In an Incresise In 
the flutter stability to maxlnun values of X for 6 of approximately 
30° cuad 1»0°, respectively, for plates conposed of two and four lamlnas. 
Thus, In designing composite laminated plates, proper fiber orientation 
Is In^ortant to Improved flutter stability 'or the plate. 

U.3.2 Effect of I.iplane Normal and Shear Loads 

The effect of the flutter boundsorles of Inplane normal and sheaur 
loads for a square, symmetric, boron-epojQr panel composed of four 
lamlnas Is shown In figures l8 and 19 where X Is presented as a func- 
tion of the Inplane load. Flutter boxindarles are shown for vedues of 
0 ® 0°, 15^. and 30°. The circle symbols Indicate the point of buckling 
with air flow. 

Inplane normal loads (See fig. l8.) result In approximately a linear 
reduction in flutter stability for loads up to the point of buckling. 

The effects of the cross-stiffness terms and rotation of the fibers away 
from the x-axls both result in a reduction In the flutter stability and 
buckling loads. Thus, classical plate theory, which neglects the cross- 
stiffness terms, is nonconsezvatlve in its predictions of both flutter 
stability and buckling loads. 

Small Inplane shear loads (See fig. 19.) have a stabilizing effect 
on flutter for G ■ 15° and G ■ 30°, but larger loads decrease the 
flutter stebillty. The Increase in flutter sterility with small inplane 
shear loads is contrary to results obtained for isotropic plates (See 
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ref. Ul.) and la due to a change in the modes that coalesce to produce 

flutter. This is shown in figure 20 where coalescence of the frequencies 

for 6 « 15° are shown for values of inplane shear loads K • 0 . , . 5 , 

xy 

and 1.0. For K >0 (See fig. 20(a).) and K ■ 1.0 (See fig. 20(c).), 
xy 

the two lowest frequencies which correspond to the (1,1) and (1,2) nx>des 

coalesce and result in flutter. However, for a smn?l. inplane shear load 

K of 0.5 (See fig, 20(b).), the lowest frequency and one of the higher 

xy 

frequencies, which correspond to the (l,l) and (2,1) modes, coalesce to 
produce flutter. For K greater than 1.0, the (1,1) and (1,2) raodt 
continue to coalesce to produce flutter but at lower values of X (See 
fig. 19.). The cross-stiffness terms and rotating the fibers away :*rom 
the x-axis both res\ilt in increased buckling loads and in improved 
flutter stability for most values of inplane shear loads. This in con- 
trary to results noted for inplane normal loeis. It should be noted 
that at the point of buckling, the flutter values of X for inplane 
shear loads (See fig. I9.) are considerably lower than those of inplane 
normal loads (See fig. I8.). This is in contrast with results presented 
in reference Ul for isotropic plates. 

4,3.3 Effect of Cross-flow 

The >‘ffects of cross- flow on the flutter boundaries for a square, 
symetric, boron-epoxy plate composed of four laminas are shown in 
figure 21. Flutter boundaries are presented for d » 0°, 15° and 30° 
where X is shown as a function of the cross-flow angle. For 6 « 15° 
and 6 * 30°, flow at small cross- flow angles may have a substantial 
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stabilizing effect on flutter. For 6 * 15°» a 30 percent increase in 
X may be obtained for a 10° cross-flow angle whereas larger cross-flow 
angles result in a reduction in flutter stability. Similar results are 
obtained for 6 = 30° where a ^^5 percent increase in X is obtained 
for a cross-flow angle of 20°. This increase in the flutter stability 
with cross-flow angle can be explained as follows: Cross-flow is de- 

stabilizing f^r an orthotropic plate (See ref. UO.). Orienting the 
fibers at an angle with the plate axis res'jlts in the maximum bending 
stiffness and thus the maximm flutter resistance occurring in that 
direction. These factors tend to counteract each other and result in 
the maximum flutter stability occurring at an angle between the x-axis 
and the fiber directions. Note that the maximuLi value of X for the 
cases shown are obtained at w wgle approximately two-thirds of the 
angle A the fibers maKe with the x-axis. 

For plates with the fibers edigned with the x-axis (6 = 0°), cross- 
flow is destabilizing, but at smsdl cross-flow angles they have higher 
values of X than plates with the fibers rotated at 0 * 15° and 6 * 
30°. Classical plate theory does not show tiny beneficial effects of 
cross- flow and thus, underpredicts the flutter stability for large cross- 
flow tmgles. For laminated con 5 >osite plates, cross-flow may be either 
stabilizing or destabilizing depending on the flow angle and the orienta- 
tion of the fibers. Thus, in designing composite plates that will experi- 
ence cross-flow, Eui optimum orientation of the fibers may considerably 
improve the flutter stability of the plate. 
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U.lt Angle-ply Plates 

Flutter boundaries calculated vising the present analysis, classical 
plate theory, and reduced bending stiffness theory are presented in this 
section for typical angle-ply plates. The effects of inplane normal and 
shear loads and cross-flow on the flutter of a square boron-epoxy plate 
composed of four laminas are also shown. There is a bending-extensional 
coupling between the governing equations for an angle-ply plate which is 
neglected by classical plate theory. Since the cross-stiffness terms 
discussed in section U.3 for symmetric plates are zero for angle-ply 
plates, the only difference between classical plate theory and the pre- 
sent ansilysis is due to bending-extensional coupling. The reduced bend- 
ing stiffness theory accounts for the coupling in an approximate way by 
reducing the bending stiffness of the plate by an amovint determined from 
the coupling terms and then neglecting the coupling in solving the equa- 
tions. For 6=0° or 6 = 90°, the coupling terms are zero, and the 
present analysis and the reduced bending stiffness theory become identi- 
cal with cleissical plate' theory. 

H.U.l Flutter Boundaries 

Flutter boundaries for square, angle-ply plates with properties 
typical of a glsuss-epoxy and a boron-epoxy material are presented in 
figure 22(a) and 22(b), respectively, where X is shown as a function 
of 0. Bending-extensional coupling has a large destabilizing effect on 
the plate composed of only two lamine^, but as the number of laminas 
increase, the coupling effect becomes smaller, and the boundaries 



approach those given by clcissicai plate theory. The reduced bending 
stiffness theory gives flutter boundaries that are in good agreement 
with the actual boundaries even for the plate composed of only two 
laminas. The small errors associated with using the reduced bending 
stiffness theory, however, are not alweys conservative. Conparing the 
flutter boundaries for the glass-epoxy (See fig. 22(a).) and the boron- 
epoxy (See fig. 22(b).) plates shows that increasing the orthotropy of 
the lamincis increases the destabilizing effect of bending-extensicneul 
coupling. The good agreement between the present analysis and the re- 
duced bending stiffness analysis does not deteriorate as the lamina 
oirthotropy is incret\sed. 

Flutter boundaries similar to those presented in figxire 22 for a 
square plate are presented in figure 23 for a boron-epoxy plate with a 
length-width ratio of 2.0. The flutter boundaries show larger destabil- 
izing effects due to bending-extensional coupling than found for the 
square plates, but increasing the number of laminas composing the plate 
reduces the destabilizing effects. For plates with four or more laminas, 
orienting the fibers at an angle wjth the x-axls may resiilt in large 
improvements in the flutter stability with the hipest flutter values 
being obtained at 6 - U5^. For plates with two laminas, rotating the 
fibers away from the x-etxis results in little or no improvement in the 
flutter stability. In each case, orienting the fibers transverse to the 
x-axis results in the most unstable condition. Good agreement is obtained 
between the present analysis and the reduced bending stiffness theory 
for all plates, but the largest discrepeuicy is obtained for the plate 



con^osed of two laminae 


U.U.2 Effect of Inplane Normal and Shear Loads 

The effects of inplane normal and shear loads on the flutter of 
square angle-ply plates are shown in figure 2U(a) and 2U(b), respectively. 
Flutter boundaries are shown for values of 0=0°, 15°, and 30°. The 
circle symbols indicate the point of buckling with air flow. Inplane 
normal and shear loads result in a sharp drop in the flutter stability 
for loads up to the point of buckling. Bending-extensional coupling has 
a destabilizing effect on flutter and reduces the panel buckling loads. 
For plates with inplane normal loads (See fig. 2k(a).), rotating the 
fibers away from the x-axis generally has a destabilizing effect on the 
flutter and reduces the buckling loads, but for plates with large in- 
plane shear loads (See fig. 2U(b).), the opposite effect occurs. How- 
ever, for small values of inplane shear, rotating the fibers away from 
the x-axis may also be destabilizing. The reduced bending stiffness 
theory gives flutter boundaries that are in good agreement v the pre- 
sent analysis for each case shown. It shoxild be noted that at the point 
of buckling, the flutter values of X for inplane shear loads are coih 
siderably lower theui those for inplane normal loads. This is in agree- 
ment with results shown for symmetric laminated plates. 

U.U.3 Effect of Cross-Flow 

The effect of cross-flow '■n the flutter boundaries for a squaire, 
angle-ply, boron-epoxy plate con^josed of four laminas is shown in figure 
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25. Flutter boundaries are presented for 6*0°, 15°, and 30° where X 
is shown as a function of the cross-flow euigle. In each case, cross- 
flow and bending-extensional coupling have a large destabilizing effect 
on the flutter boundary. However, the fibers away from the x- 

axis results in cross-flow having a less destabilizing effect ou the 
flutter boundary. Thus, for angle-ply plates, aligning the plate axis 
with the flow results in the most stable condition, but if cross-flow is 
unavoidable, proper orientation of the fiber directions may reduce its 
destabilizing effects. For each condition shown, the reduced bending 
stiffness theory shows good agreement with the present analysis. 

4.5 General Laminated Plates 

The term "general laminated plates" as iised in this study encom- 
passes plates composed of any nunber of laminas, with arbitreuy thick- 
ness, stacking sequence, and material properties, that satisfy the 
linear smedl deflection theory assumptions. This definition covers an 
infinite nuinber of plates Including the symmetric uid angle-ply plates 
discussed in sections 4.3 and 4.4. 

Flutter boundaries will be shown for two square, boron-epoxy panels 
that are composed of four laminae stacked in the sequence shown in figure 
10 and designated as P-1 and P-2. In addition, flutter bovindaries will 
be ceLLculated for six square, composite-stiffened aluminum plates con- 
figured as shown in figure 11, The flutter boundaries calculated for 
plates P- and P-2 using the present analysis will be coni)ared with 
those calculated using classiccd plate theory, and reduced bending 



stiffness theory and with those presented for the symmetric and angle- 
ply plates in sections U.3 and U.U. The difference hetween the present 
analysis and classical plate theozy for general laminated plates is due 
to both cross-stiffness terns aund bending-extensioned coupling. 

U.5.1 Flutter Boundaries 

Flutter boundaries ccdculated using the present analysis are com- 
pared with those calcvilated using classical plate theory and reduced 
bending stiffness theory in figure 26(a) and 26(b), respectively, for 
plates P-1 and P-2. The bending-extensionad coupling auid cross-stiffness 
terms have a large influence on the flutter boundary for both plates and 
render the classical plate theory hi^ly nonconservative. For these 
caises, neglecting the coupling and cross-stiffness terms say result in 
flutter values that aure close to 100 percent too high. Hovever, flutter 
boundaries calculated using the reduced bending stiffness theory show 
good agreement with the present ainadysis. The agreement is slightly 
better for plate P-1 than for P-2. This is expected since plate P-1 is 
more nearly symmetric amd thus, the coupling effect is smaller thaui for 
plate P-2. 

U.5.2 Effect of Plate Construction 

Flutter boundaries for plates P-1 and P-2 are compared in figure 27 
with those obtained for the symmetric and angle-ply plates. The curves 
shown represent flutter boundaries for six plates, all of which are com- 
posed of the same materiad amd have the same thickness but which have 
different flutter characteristics as a result of the number of layers 



and the sequence in which they are laminated. Althou^ the flutter 
hoimdaries for all the plates follow the same hasic trends with 0, the 
angle-ply plate with K U shows considerably more resistance to 
flutter them einy of the other plate constmctions. Ibe symmetric plate 
with K = 4 shows the next best resistance to flutter followed by the 
general laminated plates P-1 emd P-2 which have similar flutter boun- 
daries. The symmetric and angle-ply plates with K - 2 show the lowest 
values of X, which are as much as 40 percent lower than the maxlmxm 
vedues obtained for the angle-ply plate with K = 4. Thus, ijnprovements 
in the flutter stability are obtained by increasing the nuaiber of laminas 
and stacking them in an angle-ply sequence. 

The effects of cross-flow on the flutter boundaries for the general 
laminated plates may be seen in figure 28 where X is plotted as a fac- 
tion of the cross-flow angle. The flutter boundsuries for the symmetric 
and angle-ply plates with K * 4 are also shown for comparison. For 
the general laminated plates, cross-flow at small angles sli^tly improve 
the flutter stability. For plate P-1, cross-flow at angles greater than 
20° and for plate P-2, cross-flow at angles greater than 5° each results 
in a steady decrease in flutter stability. Similar results are obtadned 
for the symmetric plate where large Improvements in the flutter stability 
are obtained for cross-flow angles up to 10°. However, for the angle-ply 
plate, even a small amount of cross-flow results in a reduction in flutter 
stability. Although it was pointed out previously that for no cross-flow, 
the angle-ply construction is the most resistant to flutter, the presence 
of cross-flow m*^ result in the symmetric construction giving better 
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flutter stability. Thus, in designing laminated plates, the stacking 
sequence may have a significcmt effect on the flutter stability. 

U.5.3 Con 5 )osite Stiffened Alumimam Plates 

Flutter boundaries are presented in figure 29 for six square, com- 
posite stiffened aluminum plates constructed as shown in figure 11 and 
with material properties as given in table II. Two of the plates are of 
symmetric construction witn one or two layers of composite material 
applied to each side of the aluminum. The other four plates are of un- 
symmetric construction with one and two layers of conq>osite material 
applied to one or both sides of the aluminum. In each case, 50 percent 
of the plate mass is coiiq>osed of composite materials, and the only 
difference between the plates is the vay in which the composite material 
is applied. The flutter boundaries for the composite stiffened plate 
are conpared with the flutter boundary for an equal mass aluminum plate. 
The values of X shown in figure 29 are based on of the aluminm 

oi 

0 plate so that the relative advemtages of the various constructions may 

be seen. 

The composite stiffened aluminum plates have considerably hi^er 
flutter stability than the aluminum plate. For small values of 6, 
aluminim plates with laminas on both sides of the plate have values of 
. X that are 50 percent hl^er theui those obtained for plates with the 

laminas only on one side and over three times those for the eduminum 
plate. However, for increasing values of 6, the values of X decrease 
until little or no inprovement in stability is obtained by using the 
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composite material. For values of 6 between 0*^ and 90 ° ^ increasing 
the number of layers of composite material also increases the flutter 
stability. 

These results show that composite materials may be used very effec- 
tively to increase the flutter stability of an isotropic plate. For 
square plates, aligning the fibers with the plate axis and applying the 
co]]q>osite matericd symmetrically about the midplane of the plate both 
result in improved flutter stability. If the fibers are rotated with 
respect to the plate axis, increasing the number of laminas also in- 
creases the flutter stability. 
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Chapter 5 

CONCLUDING REMARKS AND RECOMMENDATIOHS 


For general laminated plates, the bending and extensional governing 
eqxiatlons are coupled and have cross -stiffness terms which are not in- 
cluded in classlcsil orthotropic plate theory. The coupling and cross- 
stiffness terms occur as a result of the lamina principal directions not 
coinciding with the plate axis. These additional, terms euid the coupling 
increase the difficulty of obtaining solutions. However, a solution 
procedure has been developed using linear small deflection theory for 
the flutter of arbitrarily laminated simply supported plates. The ex- 
tended Gederkln method is used to obtain solutions to the governing 
equations, and the aerodynamic pressure loading used in the analysis is 
that given by linear piston theory with flow at arbitrary cross-flow 
angles. 

Flutter solutions were obtained for typical symmetric, angle-ply, 
and general laminated con^oslte plates, and a limited parametric study 
was conducted. The parameters studied include the nuober, orientation, 
smd orthotropy of the lamina; the plate length-width ratio; the inplane 
normal and shear loads; end the cross-flow angle. In addition, flutter 
solutions for several conqposite stiffened aluminum plate designs were 
obtained to determine the most flutter resistant design. 

The bending-extensional coupling and the cross-stiffness terms both 
have a large destabilizing effect on the flutter of unstressed laminated 
plates, but increasing the nxuober of laminas, reducing the lamina 
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orthotropy, emd stacking the laninas in the "best" order reduce the de> 
stabilizing effect. For a squcure plate, aligning the fibers with the 
X-axis results in the hipest flutter stability, but for a plate with a 
length-width ratio of 2.0, large improvements in flutter stability may 
be obtedned by rotating the fibers away from the x-axis. For angle-ply 
plates, inplane normal and shear londs and cross- flow have a destabiliz- 
ing effect on flutter similar to that obtained for orthotropic plates. 
However, for symmetric plates with the fibers not aligned with the 
x-axis, the cross-stiffness terms result in small inplane shear loads 
and cross-flow angles inq>roving the flutter stability. Flutter calcula- 
tions for equivalent symmetric, angle-ply, and general unsymmetric 
plates indicate that for no cross- flow and no Inplane shear loads, 
plates with an angle-ply construction will have the hij^est flutter 
stability. If cross-flow or inplane shear loads are present, synasetric- 
ally constructed plates may ha\"e higher flutter stability. 

Since classical plate theory does not consider bending-extenslonal 
coupling and cross-stiffness terms, it gives Inaccurate and usually non- 
conservative flutter bovindaries for laminated plates. Reduced bending 
stiffness theory, an approximate flutter theory which accounts for the 
coupling by reducing the plate bending stiffness as determined by the 
coupling terms and then neglects the coupling in solving the equations, 
gives flutter solutions that are adeqxiate for edl plates for which numer- 
Iced results were obtained. For angle-ply plates, reduced bending stiff- 
ness theory results are obtained using published classical orthotropic 
plate theory solutions. However, for symmetric or general laminated 



plates, a solution procedure similar to the present one for symmetric 
plates is necessary to obtain even approximate solutions. 

The flutter stability of composite stiffened alumlmaa plates was 
considerably better thsm the flutter stability for an equal mass ELlumlnum 
plate. Applying the composite material to both sides of the aluminum 
plate results in better flutter stability than by applying all the 
material to one side of the plate. Also, if the fiber directions do not 
coincide with the plate axis, increasing the nusiber of layers of material 
improves the flutter stability. 

Since only a limited paraiaetric study was conducted in this investi- 
gation, it would be beneficial to use the present analysis to conduct a 
pareunet* Ic study in greater depth. Specifically, a wider variety of 
lamina materled properties should be considered and additional plates 
with general unsymmetr.ic construction should be investigated. Also, the 
possible benefits of using composite materials to stiffen conventional 
orthotropic and isotropic plates should be studied for several additional 
plates. For all results presented, the fiber directions of the laminas 
were rotated by the same absolute angle about the plate axis. Thus, 
further investigations should be made into the flutter characterl' ics 
for laminated plates in which the fib< r directions of the individual 
lamina are rotated independently. This becomes of special significance 
for laminated plates with inplane shear or cross-flow where significant 
improvements in flutter stability may be obtained by rotating the fibers 
away from the plate axis. 

The flutter results presented in the present analysis were obtained 
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for simply supported plates with free normal and tangential Inplane dis- 
placements. Other boundary conditions, which restrict the inplane dis- 
placements have been considered in the literature for simply supported 
plates and have resvOLted in large changes in natural frequencies. Since 
changes in flutter characteristics are usually associated with changes 
in the natural frequencies, the inplane boundary conditions may have a 
significant effect on the flutter characteristics. Ihus, an extension 
of the present analysis to consider other inplane boundary conditions 
would be a worthwhile endeavor. 

Although the reduced bending stiffness theory adequately predicted 
the flutter characteristics for all the plates studied, it was shown in 
the literature not to give natviral frequencies and deflections as 
accurately as desii^d for certain plates, 'nius, additional cong>arisons 
between reduced bending stiffness theory and the present analysis are 
needed for a wider range of plates. Also, since the redxiced bending 
stiffness theory does not account for the Inplane bovindary conditions, 
any change in the flutter characteristics due to the changes in boundary 
conditions suggested in the above paragraph would not be shown with the 
reduced bending ctiffness theory. These additional studies are needed 
before the reduced bending stiffness theory can be used with confidence. 

Since the literatuj*e survey revealed no nxamerical results available 
for general laminated plates, a significant contribution to the litera- 
ture could be msAe by using the present analysis to calculate the nat\aral 
vibration frequencies and the inplane normal and shear static buckling 
loads for general laminated plates. Althou^ this study has been 


I 


] 


! 


73 

concern<%d primarily with flutter characteristics « the analysis and re- 
sulting computer program can be ’ised without modifications to perform 
the inijcated calculations. 
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TABLE 

I 


LAMINA MATERIAL PROPERTIES USED IN ANALYSIS 


Material 


Properties Boron-epoxy Glass-epoxy 

30,000,000 psi 7,500,000 psi 

Egg 3,000,000 psi :2, 500, 000 psi 

1,000,000 psi 1,000,000 psi 


30 


.30 
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TABLE II 

MATERIAL PROPERTIES FOR COMPOSITE STIFFENED 
ALUMINUM HATE 


Material 


Properties 

Aluminum 

Composite 

=11 

10,000,000 psi 

29,600,000 psi 

CM 

CM 

10,000,000 psi 

U, 130, 000 psi 

®12 

3,850,000 psi 

2 , 020,000 psi 

CM 

.30 

.21 
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TABLE IV 

COMPARISON OF POUR LOWEST NATURAL FREQUENCIES 
FOR A SQUARE ANGLE-PLY PLATE ®i2^®22 * ^^12 ' 
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TABLE V 


COMPARISON OF FOUR LOWEST NATURAL FREQUENCIES 
FOR A SQUARE CROSS-PLY PLATE ^ 




\l^^P 2 

Extended Galerkin Method 

1 

6.155 

15.11*8 

15.11*8 

pl*.6l9 

10 

7.932 

21.27I+ 

21.271* 

31.727 

PO 

9.152 

25.136 

25.136 

36.608 

30 

10.210 

28.1+03 

28.1*03 

1+0.81*0 

1+0 

11 . 161 * 

31.312 

31.312 

UI4 . 656 


Rayleigh-Ritz Method, 

ref. 15 

1 

6.155 

15 . 11*8 

15.11*8 

Pl*.6l9 

10 

7.931 

21.275 

21.275 

31.721 

PO 

9.152 

25.135 

25.135 

36.601 

30 

10.210 

28.1*00 

28.1*00 

1*0.831 

1+0 

11.163 

31.308 

31.308 

1*1*. 61*6 
















Figure 1.- Geometry of angle-ply plates. 







Figure 3>” Lamina coordinate system. 
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Ftbure 7.“ Plate geometry and coordinate system. 
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K*4 

Cft) Symmetric pistes. 

Figure 11.- Nomenclature and geometry of composite stiffened aluminum 
pl&te . 








Ca) 9 = 0“. 


Figure 13.- Coalescence of frequencies for a symmetric, bo**on-epoxy 
square plate with K * It. 





Figure 13.- Concluded. 





Figure 13*“ Converged flutter bounderles for square reference plate 
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Cl) ) Boron-epoxy plate. 
Figure l6.o Concluded. 




Figure 17 •” Flutter 'boundaries for 'boron-epoxy, symmetric plate with 


a/b » 2.0. 


N “ = U 

X V n 


0 . 



103 



Figure l8.- Flutter boundaries^ for square, boron-epoxy, symmetric plate 

with Inplane normal loads. n *>9 « -0: K « U, 

y xy ’ 
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Figure 19 .- Flutter boundaries for sciuare, boroo-epoxy, symmetric plate 
wi;h inplane shear loads. 5^ ■ “ 0; K = U. 
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